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Abstract 

Extreme-scale alignment of quasar optical polarization vec- 
tors at cosmological scales (z < 2) is also characterized by 
the rotation of mean position angle x with A% w 30° per 
1 Gpc. For observing interval of z the total rotation angle 
acquires the value ~ 90°. We suggest the possible explana- 
tion of the half of this rotation as a consequence of physical 
transformation of initially vertical magnetic field By , directed 
along the normal N to the surface of accretion disk, into the 
horizontal (perpendicular to N) one. We found asymptotical 
analytical expressions for axially averaged polarization degree 
p and mean position angle x f° r various types of magnetized 
accretion disks. We found also that during the evolution can 
be realized the case B±_ w Bn where position angle x rotates 
from 45° to zero. This rotation may occur during fairly great 
cosmological time (corresponding to Az ~ 1 — 2). The part 
of rotation ~ Ax ~ 45° can be explained by a mechanism 
of alignment of polarization vectors, say distribution of the 
part of quasars as a spiral in the cosmic space with slow vari- 
ation of rotation axis of corresponding accretion disks. Both 
mechanisms are mutually related one with another. 

Keywords: polarization, magnetic fields, accretion disks, 
quasars, active galactic nuclei. 



1 Introduction 

Large scale alignments of quasar p olariz ation vectors 
have been revealed by iHutsemekerd ( 19981 ). looking at 
a sample of 170 QSO s elected from various surveys. 
Hutsemekers fc Lamvl ( 2001 ) have confirmed this effect later 
on a larger sample. The departure from random orientation 
has been found at fairly well significance level. Hutsemekers 
and Lamy have concluded that these alignments seemed to 
come from high redshift regions, implying that the underly- 
ing mecha nism might cover physical distances of gigaparsecs. 
Moreover, iHutsemekers et al. (2005) have estimated the ro- 
tation of mean position angle magnitude as Ax < 30° at the 
distance of ~ lGpc. In this paper they presented the analy- 
sis of the alignment effect for a total sample of 355 quasars, 
comprising new polarization measurements both from observ- 
ing runs 2001-2003 and new com prehensive data from s urveys 
and the literature. The paper of iBorquet et al.l (|2008l ) threw 
light on the new observational fact. They found a significant 
correlation between the polarization position angle and the 



position angle of the major axis of the host galaxy extended 
emission. 

The first question arises concerning contamination of in- 
terstellar dust grains in our Galaxy. The linear dichroism 
of aligned interstellar dust grains in our Galaxy produces re- 
ally linear polarization along the line of sight. This polariza- 
tion contaminates to some extent the quasar me asured data 
and m ay change their position angles. However. ISluse et al 



(2005) lave shown that interstellar polarization has a little 



effect on the polarization angle distribution of significantly 
polarized (p > 0.6%) quasars. 

The first explanation of observed large-scale alignment of 
quasar polarization vectors was naturally associated to light 
propa gation into intergalactic medium. IHutsemekers et ah 
(2005) claimed that, while interpretations like a global 
rotation of the Universe can potentially explain this ef- 
fect, the properties of observed phenomenon correspond, 
at first glance, to the dichroism and birefringence pre- 
dicted, for example, by photon-pseudoscalar (axion-likc 
particle) oscillation within a intergalactic magnetic field. 
This process has been considered in detail in a number 



of papers (Jain. Narain & SaralaL 20031 Das et al.. 12005 


Hutsemekers et al.l, 120051: iGnedin. Piotrovich & Natsvlishvili 


2007; IPiotrovich. Gnedin & Natsvlishvili 


^_ 2008 


Hutsemekers et al. 


. 120081 lAearwal. Kamal & Jain. 2009) 



The probability of magnetic conversion of photons into 
low- mass pseudoscalar particles was calculated in the pa- 
pers (lAnselm fc Uraltsevl. fl 982; R affelt fe Stodolskvl Il98 ' 
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lHarari fc Sikivid. I l992; Gnedin fc KrasnikovL 1 19921: iRaffel 
19961: iDeffavet et all 120011 ICsaki. Kaloper fc Terninel . |2o65 
bas et all l2005l ). 

However the basic difficulty arises because of uncertainty 
of photon-pseudoscalar (axion) mixing constant. Unfor- 
tunately, the famous CAST experiment gives only con- 
straints on this mixing ( Andriamonie et all 2007) . Recently 
IPiotrovich. Gnedin fc Natsvlishvilil (|2008i ) have estimated the 
photon-pseudoscalar mixing constant from the effect of large- 
scale alignment and rotation of polarization plane of distant 
(z > 0.5) quasars. Many authors have considered process for 
generating magnetic fields of the cosmological interest. Their 
conclusions allow us to consider unclear the magnetic field 
strength interpretation of this results, in particular, whether 
the rotation of polarization angle is associated with inter- 
galactic medium, or with the presence of some intrinsic phys- 
ical phenomena in the nearest environments of the polarized 
quasars. 

Here we present another explanation of the effect of 
extreme-scale rotation position angle versus cosmological 



1 



redshift z. We use the results of calculation of polarization 2 Basic equations 

of radiation from optically thick magnetized accretion disks 
in quasars and active galactic nuclei. The presence of 
intrinsic magnetic field in an accretion disk produces a 
new effect, provided by the Faraday rotation of polariza- 
tion plane along a photon mean free path in scatterin, 



mediu m (see, for example, iDolginov, Gnedin fc Silant 



:mg 



Silant'ev (2002) obtained the analytical approximate formu- 
lae for the Stokes parameters of polarization of radiation 
emerging from optically thick accretion disk which for the 
Milne problem acquire the form: 



(Il995h: iGnedin fc Silant'evI (Il997h : I Silant'ev! (l2002h: 



Gnedin. Silant'ev fe Shterninl (|2006h ). 

Due to these authors, a nontrivial wavelength dependence 
of polarization and position angle arises when the Faraday 
rotation angle \& at the Thomson optical depth r is sufficiently 
large: 



<£ = 0.4 
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Here, A is the radiation wavelength and 9 is the angle be- 
tween the line of sight n and magnetic field B. Below we shall 
use the known Milne problem which corresponds to the case 
where the thermal sources are located far from the surface of 
optically thick accretion disk. The existence of Faraday rota- 
tion angle ^ gives rise to both the decrease of degree of polar- 
ization p and the rotation of position angle x ( U P to maximum 
value x — 45° as compared the usual Thomson scattering 
problem where x = is assumed). The numerical solution of 
the Milne problem in the case of magnetized atmosphere with 
magnetic field directed along the normal to a t mosphere has 
been early obtained by ISilant'evj (If 994L l2002h: lAgol fc Blaesl 
<|l996h and lShternin. Gnedin fc Silant'ev! (l2003h . 

Magnetic field in an accretion disk consists of two mutually 
perpendicular components B = B 2 + B^. Here B 2 = B| is 
directed along the normal N to the disk's surface. The field 
B„ + B v is perpendicular to N, and consists of the 



Bj 

azimuthal B^ and radial (in cylindrical system of references) 
Bp components. 

We observe the azimuthally averaged Stokes parameters 
(Q(n,B)) and (C/(n,B)). If B|| = 0, the observed position 
angle x = due to symmetry of a problem (remember that 
tan 2% = U/Q, and the direction with x = is parallel to 
disk's surface). In contrary, the case Bj^ = gives rise to the 
rotation of position angle x from zero at B» = up to max- 
imum value Xmax = ±45° for large values of Bu . The signs 
± depend on outside or inside the accretion disk orientation 
of magnetic field Bp. In both cases the increase of magnetic 
field value decreases the polarization degree p. 

Note that Faraday rotation of polarization plane in inter- 
galactic medium is too low to be considered as a possible 
mechanism of observed rotation of mean position angle. 

The goal of the paper is to show that in general cases there 
exist situations where the position angle x decreases, with 
the increase of B^, from initial value x ~ 45° (corresponding 
to large Bh) up to x ~ 0. So, we shall demonstrate that the 
evolution of components B|| and Bj^ can explain the observed 
large-scale position angle rotation due to intrinsic physical 
evolution of magnetic fields in accretion disks around of QSOs 
or AGNs. Below we present this in detail. We shortly discuss 
also the origin of possible mechanisms of transformation Bn 
to B ± . 
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Q = - 



F 1 



(I- v 2 ){l + C -kfi) 



2ir Jx 1 + g (f + C - k f i) 2 + (I - q) 2 S 2 cos 2 6 ' 



U= - 



F l-g 



(f - ^i 2 )(l - q)Scost 



2tt Ji 1 + g (1 + C - k^) 2 + (1 - q) 2 S 2 cos 2 9 ' 
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where 9 is the angle between the directions of magnetic field 
B and the line of sight n, fi = cos i (i is inclination angle, 
i.e. angle between normal N and n), q is the degree of true 
absorption: q = a a /(cr a + a s ). The function J(/i) describes 
the angular distribution of emerging radiation, F is the ra- 
diation flux emerging from the disk's surface. The function 
J(l i) and the nume rical parameters g, J\, and k are tabulated 
bv ISilant'ev (2002). The minus signs denote that Thomson 
polarization is perpendicular to the plane (nN), i.e. has po- 
sition angle x = 0. In the case of dominant electron scatter- 
ing inside an accretion disk: g = 0.83255, J\ — 1.19400, and 
k = 0. Dimensio nless parameter C arises in turbulent magne- 
tized plasma (see Silant'ev ( 20051 )1 and characterizes the new 
effect - additional extinction of parameters Q and U due to 
incoherent Faraday rotations in small turbulent eddies. 

The parameters of the Faraday depolarization for Bu and 
B± can be introduced from Eq.|T|): 
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The wavelength A rest is derived in the rest system of quasar 
or AGN. The value 5 cos 9 in these notions has very simple 
form: 



S cos 9 — a + b cos $ , 



(4) 



where dimensionless parameters a and b are connected with 
the parameters 8\\ and S± 



(l-q)S\\n, b = (1 - 9 )<Wl - M 2 



(5) 



with the angle <!> = + 0*, where (f> being the azimuthal 
angle of radius-vector of observed point r on the surface of 
an accretion disk. The angle 0» is the angle between Bi and 
Bp. The azimuthal angle of line of sight n is taken zero. 

Non-polarized light escape the optically thick disk basically 
from the surface layer with tw 1. The additional extinction 
of parameters Q and U in turbulent magnetized atmosphere 
means that linearly polarized light escape the disk surface 
mainly from the Thomson optical depth t«1/(1 + C). If the 
Faraday rotation angle \& corresponding to this optical length 
becomes greater than unity, then the emerging radiation will 
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be depolarized as a result of summarizing of radiation fluxes 
with very different angles of Faraday rotation. For directions 
that are near perpendicular to the direction of the magnetic 
field in an accretion disk the Faraday rotation angle is too 
small to yield depolarization effect. Certainly, the diffusion 
of radiation in the inner parts of a disk depolarizes it even in 
the absence of magnetic field because of multiple scattering 
of photons. The Faraday rotation only increases the depo- 
larization process. It means that the polarization of outgo- 
ing radiation acquires the peak-like angular dependence with 
its maximum for the direction perpendicular to the magnetic 
field. The sharpness of the peak increases with increasing 
magnetic field strength. The basic region of allowed angles 
appears to be ~ 1/5. Another very important feature charac- 
terizing the polarized radiation is the wavelength dependence 
of polarization degree that is strongly different from that for 
Thomson scattering. 

Formulae give rise to the following expression for the 
polarization degree: 



pO*,b) 



y/(i + c- M 2 + (i - q) 2 $ 2 cos 2 e 



(6) 



where p(fx, 0) means the polarization degree for pure Thomson 
scattering. Remind that p(/x, 0) in conservative atmosphere 
(q=0) has the classical Sobolev-Chandrasekhar value, that is 
maximal and equal to 11.7% for /i = (i.e. for the incli- 
nation angle i = 90°). We use the reference system with 
X-axis lying in the plane (nN). In this system U(p,,0) = 
and p(/u, 0) = |Q(/z, 0)|. The position angle \ °f emerging 
radiation, as usually, is described by the known relation: 



U {l~q)6cos9 
tan2x= 0= 1 + C-Jfe/x- 



(7) 



For strong magnetic field strength (or large wavelength,) 
when (1 — q)5cos8 ^> 1, the simple asymptotic expressions 
take place: 



p( M ,B) 
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(8) 



It is seen from Eq. ^ that small-scale magnetic turbulence 
(parameter C) decreases the observed polarization degree. 
This is because the polarized light (parameters Q and U) 
escape mainly from the level r w 1/(1 + C) where intensity 
is lesser as compared with the level t rs 1, corresponding 
to escape of non-polarized light. In contrast, the existence 
of absorption (q ^ 0, k ^ 0) increases polarization both due 
to existence the parameter k in the denominator of Eq. (6) 
and due to that p(p, q, 0) is higher than p(/j,, q = 0, 0) (see, 
for example, Silant'ev (1980)). The absorption gives rise to 
more sharp intensity (along the normal N), and the situation 
look like the single scattering of a light beam in the surface 
layer of the atmosphere. So, at q = 0.1 the Milne problem has 
J(fi = 1) = 4.39 and p(ji = 0, 0) = 20.4%. The corresponding 
values for conservative atmosphere are 3.06 and 11.71%. 



3 Polarization degree and position 
angle of observed accretion disk 
radiation 

For an accretion disk the light depolarization depends on the 
geometry of magnetic field. Usually one observes the axially 
symmetric accretion disks as whole. In this case the observed 
integral Stokes parameters (Q(n,B)) and (U(n, B)) are de- 
scribed by the azimuthal averaged expressions. To obtain 
analytical formulae for (Q) and (U) we present expressions 
(2) for Q and U in a complex form: 
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p(h,o) 



(9) 



G + ia + z6cos3> 
Here and what follow we use the notion G = 1 + C — kp. The 
azimuthally averaging of this formula gives rise to expression: 
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(10) 



Here e = ib/(G + ia). If parameter 6^0 the ^-integration 
can be easily evalu ated by residue theorem (see, for exam- 
ple, Smirnovl ( 19641 )'). As a result, we obtain the following 
expression: 



(Q)-i(U)=± 



P (H,0) 



VG 2 + b 2 



2iGa 



(11) 



where sign plus corresponds to |e| < 1, and minus corresponds 
to |e| > 1. 

ISilant'evi (|2005l) has derived the next expression for the 
turbulent extinction parameter: 



C = 0.64(1 - q) n Xf eat (pm)(B' 2 )f B /3 



(12) 



Here, n <C 1 is the mean Thomson optical length of small 
turbulent eddies, the value B' denotes the fluctuating com- 
ponent of the magnetic field (B = Bo + B'). We omitted, 
for brevity, the subscript "0"in previous formulae. It means 
that the depolarization parameters a and b in Eq.(j3]) are de- 
termined only by the global magnetic field values. The nu- 
merical coefficient fs ~ 1 is connected with the correlation 
function of fluctuating components B' in neighboring points 
of turbulent atmosphere. 

It is convenient to introduce the relative polarization degree 
Prei — p(n,B)/p(^,0). By the usual way we obtain from 
Eq. fTTj) the following expressions: 

Prel = [G 4 + 2G 2 {a 2 + b 2 ) + (a 2 - b 2 ) 2 ] 1 ^ ' 

The relative polarization degree p re i does not depend on signs 
± in Eq. (fTTj) . The interesting property of expression (fT3f that 
it depends symmetrically on parameters a and &, i.e. we have 
p rel (n,a,b) = p re i{n,b,a). 

The position angle x is connected with the phase angle (j> of 
the complex expression inside the root in Eq.(|lip. Evidently 
we have 



tan( 



2Ga 



G 2 + b 2 



(14) 
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If right side in Eq. (fT4"|) is positive (this corresponds to sign 
plus in Eq. (fTTj) '), then the position angle x is determined by 
expression 

(U) 

tan2 X =^ = 



2Ga 



(G 2 + b 2 - 



V^G 4 + 2G 2 (a 2 + b 2 ) + (a 2 - b 2 ) 2 
For the negative 2Ga/(G 2 + b 2 — a 2 ) < we have 

tan 2\ = 
2Ga 



(15) 



v/G 4 + 2G 2 (a 2 + b 2 ) + (a 2 - b 2 ) 2 - |G 2 + 6 2 -a 2 | 



(16) 



Evidently the position angle \ is not symmetric function of 
parameters a and b. The observed parameters (U) and (Q) 
are positive for our choice of the reference frame and corre- 
sponds to the case where B is directed outside the accretion 
disk surface. If the magnetic field B is directed inside the 
accretion disk, then the position angle x changes its sign (pa- 
rameter U < 0). 

For a number of particular cases one can obtain fairly sim- 
ple analytical expressions. First of all, for the case of pure 
vertical magnetic field (parameter b = 0) we obtain directly 
from basic Eq.© the following expression: 



p re l(fj,,B\i) = 



1 



Vg 1 



For pure perpendicular magnetic field (a 
formulae give: 



tan2 X =- (17) 
0) our general 



p re ;(/z,Bj_) 



1 



VG 2 + b 2 



x = o. 



(18) 



The position angle x = is due to symmetry of problem in 
this case. For the case a ~ b, which approximately corre- 
sponds to equipartition B\\ = B± for i « 45°, the formulae 
for p and x acquire the forms: 



p re i([j,, a = b) 



tan2 X 



1 



[G 2 (G 2 - 
2Ga 



4a 2 )]V4 



(19) 



V'G 2 (G 2 + 4a 2 ) + G 2 

Remember that G = 1 + C — fc/i. It is easy check that for 
conservative atmosphere (k = 0) the existence of b — a in- 
creases the relative polarization degree p re i as compared with 
the case pure parallel magnetic field (a ^ 0, b = 0): 



p re i(jj,a,b = 0) / G 2 + , 



1/4 



p rel {^a = b) \G 2 (G 2 +Aa 2 )) <h (20) 

Below we present Tables 1-4, where p re i and x are given 
for different values of a, b and C in conservative atmosphere 
(? = 0). 

The detailed numerical calculation (particularly presented 
in tables 1-4) demonstrate that the existence of Bj^ can in- 
crease the relative degree of polarization p re i compared with 
the case 6 = 0. This effect takes place for a > 1 — 3 and 



Table 


1: 


The 


relative 


polarization 


degree 


Prei = 


p(n, B) /p(fi, 0) as a function of parameters 


a and b at G = 0. 


a\b 





2 


4 


6 


8 


10 


20 





1 


0.447 


0.242 


0.164 


0.124 


0.100 


0.050 


2 


0.447 


0.492 


0.271 


0.173 


0.128 


0.101 


0.050 


4 


0.242 


0.271 


0.352 


0.211 


0.142 


0.108 


0.050 


6 


0.164 


0.173 


0.211 


0.288 


0.178 


0.123 


0.051 


8 


0.124 


0.128 


0.142 


0.178 


0.250 


0.158 


0.054 


10 


0.100 


0.102 


0.108 


0.123 


0.158 


0.223 


0.058 


20 


0.050 


0.050 


0.051 


0.052 


0.054 


0.058 


0.149 



Table 2: The position angles x° as a function of parameters 



a and b at C 


= 0. 












a\b 


2 


4 


6 


8 


10 


20 























2 31.7 


19.0 


4.3 


1.7 


0.9 


0.6 


0.1 


4 38.0 


36.0 


20.7 


5.2 


2.3 


1.3 


0.4 


6 40.3 


39.7 


36.9 


21.3 


5.6 


2.6 


0.4 


8 41.4 


41.2 


40.3 


37.3 


21.6 


5.8 


0.7 


10 42.1 


42.0 


41.6 


40.6 


37.6 


21.8 


0.9 


20 43.6 


43.5 


43.5 


43.4 


43.3 


43.1 


22.1 



Table 


3: 


The 


relative 


polarization 


degree 


Prei = 


p(n, B) /p(fi, 0) as a function of parameters 


a and b at G = 1. 


o\6 





2 


4 


6 


8 


10 


20 





0.5 


0.354 


0.224 


0.158 


0.121 


0.098 


0.050 


2 


0.354 


0.334 


0.236 


0.165 


0.124 


0.100 


0.050 


4 


0.224 


0.236 


0.246 


0.186 


0.135 


0.106 


0.051 


6 


0.158 


0.165 


0.186 


0.203 


0.158 


0.118 


0.052 


8 


0.121 


0.124 


0.135 


0.158 


0.176 


0.140 


0.054 


10 


0.098 


0.100 


0.106 


0.118 


0.140 


0.158 


0.057 


20 


0.050 


0.050 


0.051 


0.052 


0.054 


0.057 


0.112 



Table 4: The position angles x° as a function of parameters 
a and b at G = 1. 



o\6 





2 


4 


6 


8 


10 


20 


























2 


22.5 


15.9 


6.6 


3.1 


1.8 


1.1 


0.3 


4 


31.7 


29.1 


19.0 


8.4 


4.3 


2.6 


0.6 


6 


35.8 


34.8 


30.9 


20.1 


9.2 


4.8 


0.9 


8 


38.0 


37.6 


36.0 


31.7 


20.7 


9.7 


1.3 


10 


39.3 


39.1 


38.4 


36.6 


32.2 


21.1 


1.9 


20 


42.1 


42.1 


42.0 


42.0 


41.7 


41.2 


21.8 
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b < a. The maximum polarization occurs at b = a (see tables 
1 and 3). This effect stems some "resonant" regions in an ac- 
cretion disk where the Faraday rotation from B|| is balanced 
by opposite rotation from a perpendicular magnetic field Bi. 
For b > a polarization decreases quickly with growing of pa- 
rameter b. Note also that for b = a 3> 1 the position angle 
X — > 22.5°. For this case the transition of x from w 45° to 
very small value occurs in the layer \a — b\ « 2 — 5, i.e. we ob- 
serve the rotation of position angle x from ±45° (depending 
on the magnetic field direction) to % » 0. 

Presented expressions and tables allow us to estimate the 
values of polarization degree and position angle. We shall 
use these expressions to explain the effect of the cosmological 
rot ation of the position ang le of polarized QSOs discovered 



rot ation or the position ang J 
by iHutsemekers et al. ( 2005h 



4 Scenario of intrinsic mechanism 

The scenario of our explanation is the following. First of all, 
we assume that at early time of the QSOs evolution, corre- 
sponding to z > 2 — 3, the magnetic field in accretion disk is 
directed along the normal N to the surface, and the perpen- 
dicular magnetic field Bj^ is practically absent. The increase 
of Bj^ inside the accretion disk can be considered as due to 
large-scale diffusion process, i.e. one can assume Bj^ ~ Dt. 
The time difference t is proportional to cosmological parame- 
ter z (remember that large z correspond to early time of the 
evolution of the Universe). The high values of B|| (a > 4) 
correspond to inclination angle x ~ ±45° of electric field of 
the emerging radiation, as compared to the usual Thomson 
position angle, which is parallel to accretion surface. The sign 
plus corresponds to the case where B|| is directed outside the 
surface. The minus sign corresponds to opposite field direc- 
tion. As a result of large-scale diffusion magnetic field B^ 
acquires the values B± ss B» , and the case b « a > 4 occurs. 
In small interval of growing of parameter b with \a — b\ ~ 2 — 5 
, where position angle of outgoing radiation changes from ini- 
tial value w 45° to final zero value, the angle of rotation 
can be roughly approximated by linear dependence on value 
b ~ t ~ z. Remember that usually position angle x is ob- 
served with fairly high error interval. So, evolution mech- 
anism of the increase of perpendicular magnetic field inside 
of accretion disk presents simultaneously the intrinsic mecha- 
nism of rotation of observed position angle, which is propor- 
tional to cosmological parameter z. Note that this mechanism 
is restricted by 45° rotation. 

The c osmological rotation of m ean position angle, discov- 
ered by IHutsemekers et all ( 2005 ) , corresponds to approxi- 
mately linear growing from zero at z — 2 to w 90° at z « 0. 
Our mechanism can explain only 45° rotation. Evidently 
there is an additional mechanism of rotation, connected with 
th e alignment of quas ar polarization vectors discovered also 
bv lHutsemekers ( 19981) . Clearly, this alignment is due to some 
type of anisotropy of QSOs distribution in cosmic space. For 
example, we may assume that a part of quasars distribution 
presents spiral type form. In this case the local position x = 
variates from one quasar to another. The simultaneous ac- 
tion of our intrinsic mechanism and variation of zero position 



X = for particular quasars, corresponding to z-parameter, 
can explain the observed effect of rotation of the mean posi- 
tion angle. 

The mentioned maximum angle of rotation x ~ 45° due to 
intrinsic mechanism may be not realized in reality. In this 
case main part of observed total rotation will be due to kine- 
matic mechanism of variations of the normals to accretion 
disks surfaces for particular quasars, with corresponding z - 
parameters. We only demonstrated above that the intrinsic 
mecha nism exists, and it c a n exp lain rotation < 45°. Accord- 
ing to IHutsemekers et al.l (|2005f ) the mean degree of polar- 
ization of all the sample of quasars is equal to 1.38%. They 
mentioned also that the alignment effect is more efficient for 
the low polarization quasars than for the high polarization 
ones. This means that we can consider the rotations of po- 
sition angle of quasars with the degree of polarization lesser 
then the value 1.38%. 

Now we give some examples of values of rotation for i = 
45°, 6 0° and 80°, ushrg our tables 1 and 2, and that, accord- 
ing to Chandrasekharl ( 1950h . the polarization degrees 0) 
for these inclination angles are 1.1%, 4.04%, and 6%, respec- 
tively. The accretion disk in the initial state (a = 4, b = 0) has 
degrees of polarization 0.27%, 0.98% and 1.45%, respectively, 
and the position angle x = 38°. When the growing perpen- 
dicular magnetic field B± gives rise to b = 2, accretion disk 
acquires the degrees of polarization 0.3%, 1.09%, and 1.63% 
with position angle x = 36°, i.e we have the rotation Ax = 2°. 
For b — 4 the corresponding values are 0.39%, 1.42%, 2,1% 
with position angle x = 20.7°, i.e. have rotation A% = 17.3° 
from the state (a = 4, b = 0), and the rotation A% = 15.3° 
from the state (a = 4, b = 2). Remember that values of b 
depend on z-parameter, and belong to quasars with various 
z, which are at different distances from an observer. Remind 
also that optical radiation escapes from accretion surface far 
from the centre of disk and magnetic field in this place is 
much lesser than that in the centre. The estimations give 
rise to the values B w 10 3 — 10 4 G in the central parts of a 
disk. The used value a = 4 correspond to B» ~ 16.5^t G for 
A = 0.55/im. It seems the process of diffusion of magnetic 
field from central part to the places where optical radiation 
arises can be slow. The increase of polarization degree with 
the grow of parameter b (up to b = a) is new effect which was 
explained in the end of the previous section. 

Thus, the observed effect of the cosmological rotation of 
polarization vectors of QSOs can be explained partly by evo- 
lution of a magnetic field in accretion disks. As a result of such 
evolution the topology of magnetic field in AGNs is changed 
and implies the transition, for example, from the predomi- 
nant vertical domain distribution to the predominant hori- 
zontal domain distribution. The ratio of domain sizes can be 
changed respect to the cosmological redshift z. 

The physical mechanism of magnetic inversion considered 
here can be the sa me one which was considered recently by 
Igumentsevl (J2009) for explaining the spectral transition of 
black holes binaries. The base for these phenomena can be 
the development of a magnetically arrested accretion disk at- 
tributable to the accumulation of a vertical magnetic field in 
a central part of this disk. The development and evolution 
of powerful jets provides also the evolution and inversion of 
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magnetic field components in an accretion disk. 

More interest ing mechanism has been recently suggested by 
Lvutikovl(|2009l ). Strong magnetic held can modify motion in 
the curved space-time of spinning black holes and change the 
stability conditions of circular orbits. Magnetocentrifugal jet 
launching from accretion disks around black holes is also con- 
nected to_top_ok)gy_of accretion disk magnetic held. Accord- 
ing to Lvutikov ( 20091 ). magnetocentrifugal launching for a 
Schwarzschild black hole requires that the poloidal component 
of magnetic held makes an angle less that 60° to the outward 
direction at the disk surface. For the prograde rotating disks 
around Kerr black holes this angle increases and becomes 90° 
for footpoints anchored to the inner region of an accretion 
disk for a limit spinning a* = 1 black hole (a* = cl/Mbh)- It 
means that the effect of cosmological rotation of QSO polar- 
ization plane can be interpreted as a result of spin evolution 
of QSO: a* w 1 for z > 1 and a* < 1 for z < 1. This 
suggestion has some evidence. Recently, the constraints on 
the spins of the black holes in the nearby QSOs have been 
obtained: a* = 0.6 ± 0.2 in the n arrow-line Seyfert SWIPT 

" 20091) o, = 0.60 ± 0.07 for the 



12127+5654 ( Miniuttietal 



SMBH in Fairall 9 (|Schmoll et all . [2009) 



The another idea has been developed bv lContopoulos et al 
(|2009h . They developed the scenario in which cosmic mag- 



netic fields are generated near the inner edges of accretion 
disks in AGNs by azimuthal electric currents due to differ- 
ence between the plasma electron and ion velocities that 
arises when the electrons are retarded by interactions with 
X-ray photons. This mechanism namely relates the polarity 
of poloidal magnetic held to the angular velocity of the ac- 
cretion disk, i.e. this polarity depends strongly on spin of a 
black hole. If the spin decreases at low redshifts z < 1, this 
effect produces magnetic inversion in the accretion disks. 

Unfortunately, if the effect of rotation of polarization planes 
in QSOs can be partly explained in framework of our mecha- 
nism, the origin of cosmological alignment of polarization vec- 
tors becomes open. It is not except, that the mechanism of 
cosmological differential rotation produces such kind of align- 
ment. It seems more probable supposition of alignment mech- 
anism is the existence of turbulent motions of the huge-scale. 

Recently, the int eresting idea of alignment of gala xies has 
been suggested by ITruiillo, Carretero fc Patr 3 (|2007h . They 
claimed that galaxies are not distributed randomly through- 
out space but are arranged in "cosmic web" of filaments and 
walls. Unfortunately, there is still no compelling observa- 
tional evidence of a link between the structure of the cosmic 
web and how galaxies form within it. The basis for this con- 
nection is the origin of galaxy angular momentum. The spin 
of galaxies can be generated by tidal torques operating in 
the early Universe on primordial material from which galaxy 
is formed. In its order, magnetic helds in QSOs and galax- 
ies can corr elate with galaxy rotation rate an d its angular 



momentum. ITruiillo. Carretero fc Patril (|2007l ) claimed that 
observational link between large-scale structure and the prop- 
erties of individual galaxies has been definitively estab lished. 
One should remind the result bv lBorquet et al. I (j2008l) . They 
found a correlation between the polarization position angle 
and the position angle of the major axis of the host galaxy 
extended emission. It seems that in the frame of this idea 



our intrinsic mechanism can explain the total rotation of the 
mean position angle if we shall suppose that position of zero's 
X in neighboring "webs" has a jump w 45°. 

5 Conclusions 



We showed that the observed by lHutsemekers et al. ( 2005h ef- 
fect of cosmological rotation of polarization vectors of QSOs 
can be partly explained by the evolution of magnetic held 
in accretion disks around SMBHs. The evolution of topol- 
ogy of magnetic held in AGNs produces the transition from 
the predominant vertical domain distribution at z > 1 to the 
predominant horizontal domain distribution at z < 1. Such 
kind of transition can be explained by the cosmological evo- 
lution of spin of supermassive black holes from the Kerr case 
at z > 1 to Schwarzschild case at z < 1. The calculations 
show that this effect is able to provide the observed rotation 
measure at the level of 30° / Gpc. The alignment effect itself 
is leaving unexplained. One of possible scenarios is alignment 
of spins of galaxies with observed large-scale structure of the 
Universe, which possibly represent the turbulent eddies of the 
huge scale. 

Acknowledgements 

This research was supported by the RFBR (project No. 07- 
02-00535a), the program of Prezidium of RAS "Origin and 
Evolution of Stars and Galaxies" , the program of the Depart- 
ment of Physical Sciences of RAS "Extended Objects in the 
Universe" and by the Grant from President of the Russian 
Federation "The Basic Scientific Schools" NS-6110.2008.2. 
M.Yu. Piotrovich acknowledges the Council of Grants of 
the President of the Russian Federation for Young Scientists, 
grant No. 4101.2008.2. 

References 

Agarwal N., Kamal A., Jain P., 2009, arXiv: 0911.0429 

Agol E., Blaes O., 1996, MNRAS, 282, 965 

Andriamonje S., Aune S., Autiero D. et al., 2007, JCAP 0704, 
010 

Anselm A.A., Uraltsev N.G., 1982, Phys.Lett.B, 116, 161 

Borguet B., Hutsemekers D., Letawe G. et al., 2008, A&A, 
478, 321 

C arilli C.L., Miley G.K ., Rottgering H. et al., 2000, 
|arXiv:astro-ph/0008380| 

Chandrasekhar S., 1950, Radiative transfer, Clarendon Press, 
Oxford 

Contopoulos I., Christodoulou D.M., Kazanas D., Gaburda 
D.C., 2009, ApJ, 702, L148 

Csaki C, Kaloper N., Terning J., 2002, Phys.Rev.Lett., 88, 
161302 



6 



Das S., Jain P., Ralston J.P., Saha R., 2005, JCAP, 06, 002 

Deffa yet C, Harari P. , Uzan J.-P., Zaldarriaga M., 2001, 
|arXiv:hep~ph/0112118| 

Dolginov A.Z., Gncdin Yu.N., Silant'ev N.A., 1995, Propaga- 
tion and polarization of radiation in cosmic media, Gordon 
& Breach, New York 

Gnedin Yu.N., Krasnikov S.V., 1992, Sov.Phys.JETP, 75, 933 

Gnedin Yu.N., Silant'ev N.A., 1997, Basic Mechanisms of 
Light Polarization in Cosmic Media, Hartwood, Amster- 
dam 

Gnedin Yu.N., Silant'ev N.A., Shternin P.S., 2006, As- 
tron.Lett., 32, 39 

Gnedin Yu.N., Piotrovich M.Yu, Natsvlishvili T.M., 2007, 
MNRAS, 374, 276 

Harari D., Sikivie P., 1992, Phys.Lett.B, 289, 67 

Hutsemekers D., 1998, A&A, 332, 410 

Hutsemekers D., Lamy H., 2001, A&A, 367, 381 

Hutsemekers D., Cabanac R., Lamy H., Sluse D., 2005, A&A, 
441, 915 

Hutsemekers D., Payez A., Cabanac R. et al., 2008, 
larXiv:0809.3088l 

Impey CD., Lawrence C.R., Tapia S., 1991, ApJ, 375, 46 

Jain P., Narain G., Sarala S., 2003, |arXiv:astro-ph/0301530| 

Igumentsev I.V., 2009, ApJ, 702, L72 

Lyutikov M., 2009, MNRAS, 396, 1545 

Miniutti G., Panessa F., de Rosa A. et al., 2009, MNRAS, 
398, 255 

Raffelt G.G., Stodolsky L., 1988, Phys.Rev.D, 37, 1237 

Raffclt G.G., 1996, Stars as Laboratories for Fundamental 
Physics, The University of Chicago Press, Chicago 

Piotrovich M.Yu., Gnedin Yu.N., Natsvlishvili T.M., 2008, 
larXiv:080O649l 

Schmoll S., Miller J.M., Volonteri M. et al., 2009, ApJ, 703, 
2171 

Sluse D., Hutsemekers D., Lamy H., Cabanac R., Quintana 
H, 2005, A&A, 433, 757 

Shternin P.S., Gncdin Yu.N., Silant'ev N.A., 2003, As- 
trofizika, 46, 433 

Silant'ev N.A., 1980, SvA., 24, 195 

Silant'ev N.A., 1994, JQRST, 52, 207 

Silant'ev N.A., 2002, A&A, 383, 326 

Silant'ev N.A., 2005, A&A, 433, 1117 



Smirnov V.I., 1964, Course of Higher mathematics, Vol 3, 
part 2, Pergamon Press, London 

Trujillo I., Carretero C, Patri S.G., 2007, Cosmic Frontiers 
ASP Conference, 379, 350 

Velikhov E.P., 1959, Sov.Phys.JETP Lett., 36, 1398 

Wills B.J., Wills D., Breger M. et al., 1992, ApJ, 398, 454 



7 



